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We discuss the new class of static axially symmetric black hole solutions obtained recently in 
Einstein- Yang-Mills and Einstein- Yang-Mills-dilaton theory. These black hole solutions are asymp- 
totically flat and they possess a regular event horizon. The event horizon is almost spherically 
symmetric with a slight elongation along the symmetry axis. The energy density of the matter 
fields is angle-dependent at the horizon. The static axially symmetric black hole solutions satisfy a 
simple relation between mass, dilaton charge, entropy and temperature. The black hole solutions 
, are characterized by two integers, the winding number n and the node number k of the purely 

■ magnetic gauge field. With increasing node number the magnetically neutral black hole solutions 

^\ ' form sequences tending to limiting solutions with magnetic charge n, corresponding to Einstein- 

Maxwell-dilaton black hole solutions for finite dilaton coupling constant and to Reissner-Nordstr0m 
' black hole solutions for vanishing dilaton coupling constant. 

Q 

(N 

> 
00 

o 

(N 

7' 

u 



X 



Preprint br-qc/9712086 



1 



I. INTRODUCTION 



The "no hair" conjecture for black holes states, that black holes are completely characterized by their mass M, 
their charge Q and their angular momentum J. This conjecture presents a generalization of rigorous results obtained 
for scalar fields coupled to gravity ^ as well as for Einstein-Maxwell (EM) theory In EM theory, the unique 
family of stationary Kerr-Newman black holes with nontrivial values of M, Q, and J contains the stationary Kerr 
black holes for Q — the static Reissner-Nordstr0m black holes for J = and the static Schwarzschild black holes 
for J = Q = Q. Notably, the static black hole solutions in EM theory are spherically symmetric, and the stationary 
black hole solutions are axially symmetric. 

In recent years counterexamples to the "no hair" conjecture were established in various theories with non-abelian 
fields, including Einstein- Yang-Mills (EYM) theory, Einstein- Yang-Mills-dilaton (EYMD) theory, Einstein- Yang-Mills- 
Higgs (EYMH) theory, and Einstein-Skyrme (ES) theory [^. Possessing non-trivial matter fields outside their regular 
event horizon, these non-abelian black hole solutions are no more completely determined by their global charges. And 
they hold more surprises. There are static black hole solutions with only axial symmetry Q|, static black hole solutions 
with only discrete symmetries and there are non-static non-rotational black hole solutions 1^. Unlike the static 
spherically symmetric and axially symmetric static black hole solutions |^ , many of the new types of non-abelian 
black hole solutions are only perturbative solutions |^|-@|. 

The static axially symmetric black hole solutions in EYM and EYMD theory Q have many properties in common 
with the globally regular static axially symmetric solutions constructed previously ,^ . Representing generalizations 
of the static spherically symmetric globally regular and black hole solutions |p^|l2|, these static axially symmetric 
solutions are characterized by two integers. These are the node number k of the gauge field functions and the winding 
number n with respect to the azimuthal angle </>. While (j) covers the full trigonometric circle once, the fields wind n 
times around. The static spherically symmetric solutions have winding number n = 1. Winding number n > 1 leads 
to axially symmetric solutions. 

The static axially symmetric EYM and EYMD black hole solutions Q are asymptotically flat and possess a regular 
event horizon. The event horizon of these solutions resides at a surface of constant isotropic radial coordinate. The 
energy density of the matter fields is not constant at the horizon but angle-dependent. Outside their regular event 
horizon, the static axially symmetric black hole solutions possess non-trivial magnetic gauge field configurations, but 
they carry no global magnetic charge. For fixed winding number n and increasing node number k the solutions form 
sequences, tending to limiting solutions. These limiting solutions are spherically symmetric and abelian, representing 
Einstein-Maxwell-dilaton (EMD) and Reissner-Nordstr0m (RN) black hole solutions for EYMD and EYM theory, 
respectively which carry magnetic charge n. 

Having given a brief account of the properties of these static axially symmetric black hole solutions in Q , we here 
present these solutions in detail. The paper thus represents the second paper of our sequel on static axially symmetric 
EYMD solutions, following [^, where a detailed presentation of the globally regular solutions was given. In section 
II of this paper we recall the action, we present the static axially symmetric ansatz in isotropic coordinates and 
we discuss the boundary conditions. Introducing temperature and entropy of the black hole solutions, we derive a 
relation between these thermodynamic quantities and the mass and the dilaton charge. In section III we recall the 
static spherically symmetric black hole solutions and present them in isotropic coordinates. In section IV we present 
the static axially symmetric black hole solutions, discuss their properties, and analyse the properties of their event 
horizon. We show that the sequences of neutral non-abelian solutions tend to limiting charged abelian solutions. We 
present our conclusions in section V. In Appendix A we present the expansion of the functions at the horizon. We 
show that the surface gravity is constant at the horizon and that the Kretschmann scalar is finite there. In Appendix 
B we discuss the final choice of functions for the numerical integration. 

II. STATIC AXIALLY SYMMETRIC ANSATZ 



We consider the SU(2) Einstein- Yang-Mills-dilaton action 



S= I (T^ + LM]V^9d'x (1) 



with the matter Lagrangian 



Lm = -^d^^d''^ - e2«*iTr(F^,F^'') , (2) 
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the field strength tensor 



Ff,^ = df,A^ " d^A^ + ie [Af,,A^] , (3) 

the gauge field 

A, - ^tM^ , (4) 

the dilaton field <I>, and the Yang-Mills and dilaton coupling constants e and k, respectively. 

Variation of the action (|l|) with respect to the metric g'^'^ leads to the Einstein equations, variation with respect to 
the gauge field and the dilaton field $ leads to the matter field equations ||] . 



A. Static axially symmetric ansatz 

As for the globally regular static axially symmetric solutions HJ^, we adopt isotropic coordinates to construct the 
static axially symmetric black hole solutions. In terms of the coordinates r, 6 and (j) the isotropic metric reads 

ds^ = -fdr + -jdr^ + —dd + J d4f , (5) 

where the metric functions /, m and / are only functions of the coordinates r and 9. Regularity on the z-axis 
{z ~ r cos 9) requires 

m|e=o = '|e=o • (6) 
Again, we consider a purely magnetic gauge field, Aq = 0, and choose for the gauge field the ansatz jlq-O, 



A^^dx^ = ^ [r^' [Hidr + (l - H2) rd9) - n [tI^H:, + (1 - Ha)) r sin 9d^] , (7) 

where the symbols r", and denote the dot products of the cartesian vector of Pauli matrices, f = {Tx,Ty,Tz), 
with the spatial unit vectors 

e " = (sin 6 cos ncj), sin 6 sin n(f>, cos 6) , 
eg = (cos 9 cos n(f>, cos 9 sin ncf), — sin 9) , 

= (— sin n0, cos ncj), 0) , (8) 

respectively. Since the fields wind n times around, while the azimuthal angle (p covers the full trigonometric circle 
once, we refer to the integer n as the winding number of the solutions. The four gauge field functions Hi and the 
dilaton function $ depend only on the coordinates r and 9. The spherically symmetric ansatz p2| is recovered for 
n = 1 and Hi ^ H3 = 0, H2 ^ H4 ^ w{r) and <i> = $(r). 

The ansatz (^ is axially symmetric in the sense, that a rotation around the z-axis can be compensated by a gauge 
rotation. Besides being axially symmetric the ansatz respects the discrete mirror symmetry M^z 'Si C, where the first 
factor represents reflection through the xz-plane and the second factor denotes charge conjugation |l^ , [l^ , [l9| . 

The ansatz is form-invariant under the abelian gauge transformation [p^ , ^^^ 

U = exp(^^T^r{r,9)^ . (9) 

The functions Hi and H2 transform inhomogcneously under this gauge transformation, 

Hi^ Hi~ rdrT , 

H2~^H2+ deV , (10) 

like a 2-dimensional gauge field. The functions Hs and H4 combine to form a scalar doublet, {H^ + ctg9,H4). We 
choose the same gauge condition as previously ||l^ , [l(j , [l7| ,p|,[i|Js|l , 

rdrHi - deH2 = . (11) 
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With the ansatz (||)-(0) and the gauge condition (|T^) we obtain the set of EYMD field equations, given in 
The energy density of the matter fields e = —Tq = —Lm reads 



-in — 



J_ 

2m 

2 



^ {rdrH2 + deH,f 



I 



ideH3 - 1 + ctgOHs + H2Hif + {deHi + ctgO {Hi - H2) - iJzFs)"] } 



(12) 



Here the first gauge field term derives from Fre, the second and third derive from Frr/, and the fourth and fifth from 
i<6(0- As seen from eq. (12), regularity on the z-axis requires 



2\e=o 



(13) 



B. Boundary conditions 



To obtain asymptotically flat solutions with a regular event horizon and with the proper symmetries, we must 
impose the appropriate boundary conditions, the boundaries being the horizon and radial infinity, the z-axis and, 
because of parity reflection symmetry, the p-axis. The presence of a regular event horizon is the essential new feature 
of the static axially symmetric black hole solutions, as compared to the globally regular solutions. We therefore 
begin with a detailed discussion of the boundary conditions at the horizon. The boundary conditions at infinity and 
along the p- and the z-axis agree with those of the globally regular solutions . They are only briefly recalled for 
completeness. 

Boundary conditions at the horizon 

The event horizon of the static black hole solutions is characterized by = — / = 0, grr is finite at the horizon 
in isotropic coordinates. We impose that the horizon of the black hole solutions resides at a surface of constant r, 
r = th [ pO| . This ansatz for the event horizon is justified a posteriori, since it leads to consistent solutions, possessing 
a regular event horizon. 

Requiring the horizon to be regular, we obtain the boundary conditions at the horizon ?' = rpj. The metric functions 
must satisfy 

/|r=rH = 'm\r-=rH = 1\t=th = , (14) 

and the dilaton function 

9.$|.=.H = . (15) 

The conditions for the gauge field functions are 

{dsHi+rdrH2)\r=r^^0 , 

{rdrHs - HiHi)\r=r^ = 0, {rdrHi + Hi{H3 + Ctg0)) \r=r^ = , (16) 

which imply, Fre = and Fr^ — 0, respectively. 

Thus the equations of motion yield only three boundary conditions for the four gauge field functions Hf, one gauge 
field boundary condition is left indeterminate. However, for the black hole solutions precisely one free boundary 
condition at the horizon is necessary to completely fix the gauge. The reason is, that in contrast to the case of 
the globally regular solutions ||], for the black hole solutions the gauge condition (|ll|) still allows non-trivial gauge 
transformations satisfying 

r'^d^r + rdrV + d'^T = . (17) 
To fix the gauge, we have implemented various gauge conditions at the horizon, such as 

idrHi)\r^rn = , (18) 



or 
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{deHi)\r^rn = , (19) 

obtaining the same results for the gauge invariant quantities. 

The expansions of the functions at the regular horizon are given in Appendix A. A. 
Boundary conditions at infinity 

At infinity (r = oo) we require the boundary conditions |^:^,^ 

f\r—oo — ^|r— oo — ^|r— oo — 1 ; (^0) 
$|r=oo=0, (21) 

(since any finite value of the dilaton field at infinity can be transformed to zero via <i) $ — $(00), r re^'^*'^°°-') 
and 

^^2|r=oo — ^^4|r=oo ~ ±1, ^^l|r=oo = ^^3|r=oo = , (22) 

to obtain magnetically neutral solutions. The expansion of the functions at infinity is given in 

The node number k is defined by the number of nodes of the gauge field functions H2 and H4 Q . Because of the 

symmetry with respect to the transformation Hi — > —Hi, we can choose these gauge field functions to be positive at 

the horizon. Solutions with an even number of nodes then have i?2(oo) — Hi^oo) — 1, whereas solutions with an odd 

number of nodes have H2{oo) — _ff4(oo) = —1. 
Boundary conditions along the axes 

The symmetries determine the boundary conditions along the p-axis and the z-axis 



(23) 
(24) 



def\e=a = 90m|e=o = del\e=o = 

def\e=n/2 = dgm\g=^/2 = dgl\g=^/2 = 

5e$|e=o = , 
de^g=^/2 = , 

and 

Hi\g=o — H^le^o = , dgH2\g=Q = dgH4\g=o — , 
Hl\g=Tr/2 — H^\e^.^/2 — , deH2\g=T:/2 = dgHi\g^j^/2 — . 

In addition, regularity on the z-axis requires condition (^ for the metric functions to be satisfied (see Appendix B) 
and condition ( |l3|) for the gauge field functions. The expansion of the functions on the positive z-axis is given in 

C. Dimensionless Quantities 

As previously, we introduce the dimensionless coordinate x, 

e 

the dimensionless dilaton function ip, 



X = -^=r , (26) 



(f = V4^$ , (27) 
and the dimensionless dilaton coupling constant 7, 

The dilaton coupling constant 7 represents a parameter; for 7=1 contact with the low energy effective action of 
string theory is made, whereas in the limit 7 — > the dilaton decouples and EYM theory is obtained. 
The dimensionless mass /i is related to the mass M via 

eG , , 
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D. Mass, temperature and entropy 



The mass M of the black hole solutions can be obtained directly from the total energy-momentum "tensor" r^'' of 
matter and gravitation pi| 



M 



(30) 



As shown in |9[] , the dimcnsionless mass /i is then determined by the derivative of the metric function / at infinity 



/i = i hm x'^dxf 

Z X — >oo 



(31) 



The zeroth law of black hole physics states, that the surface gravity n^g is constant at the horizon of the black hole 
solutions, here lE2,pl 



< = -(l/4)gV(9.ff«)(a,3«) . 



(32) 



To evaluate Kgg, we need to consider the metric functions at the horizon. Expanding the equations in the vicinity of 
the horizon in powers of the dimcnsionless coordinate {x — x-ii)/x-ii, we observe, that the metric functions arc quadratic 
in a; — .th. 



f{x,e)^h{0) 



X ^ xn 
xn 



O 



m(x, 9) = m2{9) 



X ~ xn 
xn 



O 



x-a 



(33) 



and likewise for l{x, 9) (see Appendix A. A). Since the temperature T is proportional to the surface gravity Kgg p^ , 

T = «,g/(2^) , (34) 

we obtain for the dimensionless temperature 

f2{0) 



T = 



2TTX]i\/m2{9) 



(35) 



With help of the Einstein equations we show in Appendix A.B, that the temperature T, as given in (|35|), is indeed 
constant. 

The dimensionless area A of the event horizon of the black hole solutions is given by 



A = 27r / d9 sin 



The entropy S is proportional to the area A [E2[ 



4 ' 



leading to the dimensionless product 



TS 



■EH 

4 7o 



d9 sin ( 



(36) 



(37) 



(38) 



Having defined temperature and entropy, we now derive a second expression for the mass |22|] . As shown before , 
the equations of motion yield 



SttG 



Integrating both sides over r, and from the horizon to infinity, we obtain 

^27r /"TT 1*00 



1 

AG 



d9 sin 9 







/ / / d(pd9dr^ (2To° - T/) = 

Jo Jo Jrn 



(39) 



(40) 
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Changing to dimensionless coordinates, we express the l.h.s. with help of the dimensionless mass /i and the product 
of temperature and entropy, TS, obtaining 

fi = fio + 2TS , (41) 



with fio ~ [e/ \/ A'kG)GMo, in agreement with the general mass formula for static black hole solutions [ p2[ . 
Considering the dilaton field equation in the form |^ 

-^d^ (7=^9''$) = {V^d^^ In /) , (42) 

it is straightforward to derive a relation between mass the /i and the dilaton charge D, 

D = lim x'^d^ip . (43) 



We replace the l.h.s. of eq. ( |39[) by the l.h.s. of eq. (|42| ) and integrate as above. Recalling the boundary condition 
dr^\r=rii — and changing to dimensionless coordinates, we obtain 

D = 7Mo = 7(/^ - 2TS) , (44) 

where the second equality follows from eq. (^l|). 

Previously we obtained for static spherically symmetric black holes the relation 

^ = 7(M-gr) , (45) 

where fis = 5:h/2 and Tg — l/(47ra;H) represent the Schwarzschild mass and the Schwarzschild temperature of a black 
hole with horizon xn in Schwarzschild-like coordinates. This relation represents a special case of the more general 



relation (44), since for a spherical black hole A = Ajtx^. 



III. SPHERICALLY SYMMETRIC SOLUTIONS 

The static spherically symmetric EYM and EYMD solutions were obtained previously in Schwarzschild-like coor- 
dinates with metric ||l^,|li,|3) 

ds^ = -A^Ndt^ + j^df-^ + f2 (d6|2 + sin^ dd(l>^) (46) 



and metric functions A{f) and N{r), 



Nif) ^ 1 _ 2m(f) _ ^^^^ 



We here briefly present the spherical solutions in isotropic coordinates. We exhibit the coordinate transformation 
between the radial coordinates r and f, and we discuss the limiting solutions. 

A. Coordinate transformation 

The static spherically symmetric isotropic metric reads 

ds^ = -fdf + J {dr^ + {dO^ + sin^ Odcj)^)) . (48) 

The coordinates r and r are related by 

(49) 
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Changing to dimensionless coordinates x and x, the coordinate function x{x) must be obtained numerically from 
eq. (^9|), since the function N(x) of the non-abelian solutions is only known numerically. To avoid the divergence of 
1/N{x) at the horizon in the numerical integration, we introduce the function A(a:) 



A{x) 



where 



and Ng is the Schwarzschild metric function 



.dN 
dx 



(50) 



(51) 



X 



(52) 



After replacing on the r.h.s. of cq. ( [49| ) l/\/N by 1 / c\/ Ns + A, we integrate the first term analytically and the finite 
A-term numerically. This yields the coordinate transformation 



X = 



2y^i(i — xu) +2x — Xu 



where 



exp 



(53) 



Xh 



XH 

4 



(54) 



is determined by the asymptotic requirement, x/x — > 1. For the Schwarzschild solution c = 1, and xn — xn/A. 

Fig. 1 demonstrates the coordinate transformation for the static spherically symmetric EYMD solutions with 
xh — 0.02, 7 = 1 and k — 1 — 4. Figs. 2a, b show the metric functions / and m, and Figs. 3a, b show the gauge field 
function w and the dilaton function if. 

Fig. 4 demonstrates the coordinate transformation for the static spherically symmetric EYM solutions with a:H = 
0.02 and k — 1 — A. Figs. 5a,b show the metric functions / and m. 



B. Limiting solutions 

For fixed dilaton coupling constant and horizon, the sequences of neutral static spherically symmetric EYMD black 
hole solutions converge to limiting solutions . These limiting solutions are EMD black hole solutions |jl^ with the 
same dilaton coupling constant, the same horizon and charge P — I | p3[ |. 

We now consider the limiting EMD solution for 7 = 1 and a general charge P. Defining |n3| 



2P' 



X+ = ^il + 2P2 , X_ = — , (55) 

and 



X- + ^4x2 + Xl 
2 

i.e. X-n = the coordinate transformation for the limiting EMD solution reads 



X = '-^ , (56) 



with 



2^(X - X^){X - X+) + 2X- (X_ + X+) 
X = — , (57) 



= ^i-^ . (58) 
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The metric functions of the Hmiting solution read in isotropic coordinates 



1 + 2 



/ ^ (59) 



and 

m^=(l-f^Y)\ (60) 



The dilaton function of the limiting solution reads 



(1 + my 



l + 2SiA/l + ^^) +(^ 

The gauge field function of the limiting solution is trivial, Woo = 0. For P ~ 1 and xh = 0.02 the coordinate 
transformation for the limiting solution is shown in Fig. 1, the metric functions of the limiting solution are shown in 
Figs. 2 and the dilaton function in Fig. 3. For the horizon xh = 0.02 the convergence is rapid. The solution with 
k — S is already very close to the limiting solution, and the solution with A: = 4 is almost indistinguishable, for all 
functions except for the gauge field function w. This function approaches its limiting function Woo — non-uniformly, 
since the boundary conditions require w{oo) ^ 0. A detailed discussion on the convergence of these solutions is given 
in||]. 

We now turn to the limiting solutions of the the sequences of static spherically symmetric EYM black hole solutions. 
In Schwarzschild-like coordinates, the limiting solution of the sequence with fixed horizon in, is a RN solution with 
the same horizon and with charge P = 1, if in > 1 ||2^ , p3|| . If xh < 1, the limiting solution consists of two parts, an 
exterior part covering the interval 1 < i < oo, which represents the exterior of an extremal RN solution with mass 
/i = 1, horizon = 1 and charge P = 1, and an oscillating interior part covering the interval in < i < 1 ||2^ , ^ , p5| . 
In isotropic coordinates with fixed horizon xjj , the limiting solution corresponds to the exterior of a RN solution with 
the same horizon and with charge P = 1. 

The coordinate transformation for the limiting RN solution reads 



a/x^ - 2^x + P2 + 5 - ^ 
X = , (62) 



with 



The metric functions are given by 



(63) 



/oo = ; TT (64) 



and 



-(f)) ■ 

i.e. rrioo is identical for the EMD and RN solutions. Again, the gauge field function of the limiting solution is trivial, 
Woo = 0. For P = 1 and xh = 0.02 the coordinate transformation of the limiting solution is shown in Fig. 4 and the 
metric functions in Figs. 5. 
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C. Limit xm — > 



Let us now consider the limit xh for the black hole solutions. In this limit the solutions tend towards the 
corresponding globally regular solutions ||ll|,|l^,|3| . For several quantities of interest, however, the limit a;H — + is not 
smooth. For instance the energy density of the matter fields of the black hole solutions approaches the energy density 
of the globally regular solutions with a discontinuity at the origin. The reason is, that the magnetic field of the black 
hole solutions is purely radial at the horizon 

B = BrCr , (66) 

with Br — Fg^, because the boundary conditions (|l^) require 

Bg^O , B^^Q , (67) 

with Bg = — and B^ = Frg, whereas the magnetic field of the globally regular solutions has non- vanishing Bg at 
the origin. We demonstrate this discontinuous behaviour of the energy density of the matter fields in Figs. 6a,b for 
the EYMD (7 = 1) solutions with fc = 1 and 2 and xu = 0.001, 0.01, 0.1 and 1. 
We finally consider the Kretschmann scalar K, 

K = R'^^'^^R^^o.^ . (68) 

An analytical expression for K is given in Appendix A.C. We show K in Fig. 7 for the same set of EYMD solutions. 
Again, the limit .th ^ is not smooth. However, as required for a regular horizon, K is finite at the horizon for finite 



IV. AXIALLY SYMMETRIC SOLUTIONS 



Subject to the above boundary conditions, we solve the equations for the static axially symmetric black hole 
solutions numerically. We employ the same numerical algorithm |26| as for the static axially symmetric globally 
regular solutions |^,^,^. To map spatial infinity to the finite value x = 1, we here employ the radial coordinate 

X = 1 - — . (69) 

X 

The equations are then discretized on a non-equidistant grid in x and 6, where typical grids used have sizes 150 x 30, 
covering the integration region 0<a;<l, O<0< 7r/2. The numerical error for the functions is estimated to be on 
the order of 10~^. 

The solutions depend on two continuous parameters, the "isotropic radius" xn of the horizon and the dilaton 
coupling constant 7, as well as on two integers, the winding number n and the node number k. 



A. Energy density and horizon 

We begin our discussion of the static axially symmetric black hole solutions by considering the energy density of 
the matter fields. As an example we show in Figs. 8 the energy density of the matter fields for the black hole solution 
with a;H = l,7 = l, ^i = 2 and k ~ 1. Fig. 8a shows a 3-dimensional plot of the energy density as a function of the 
coordinates p = xsinO and z — xcosO together with a contour plot, and Figs. 8b-e show surfaces of constant energy 
density. For small values of e the energy density appears ellipsoidal, being fiatter at the poles than in the equatorial 
plane. With increasing values of e a toruslike shape appears with two additional ellipsoids covering the poles. The 
ellipsoids covering the poles persist up to the largest values of the energy density, showing that the maximum of the 
energy density resides at the poles. Furthermore these black hole solutions have the remarkable property, that the 
energy density is not constant at the horizon but angle-dependent. 

The static axially symmetric black hole solutions are self-consistent solutions arising from the interplay of gravity 
with the non-abelian gauge fields. In isotropic coordinates the horizon of the static axially symmetric black hole 
solutions resides at a surface of constant radial coordinate x, x — xn. Since the energy density of the matter fields of 
the static axially symmetric black hole solutions is angle-dependent at the horizon, this suggests that the horizon is 
deformed. We therefore measure the circumference of the horizon along the equator, Lg, 
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:a;sinf 



27ra;H 



and the circumference of the horizon along the poles, Lp, 



(70) 



L„ = 2 / d6 J — 



= 2a;H / de 



(71) 



x—xh,<P— const. 

A spherical horizon would require = Lp. For the static axially symmetric black hole solutions we observe, however. 



Lp > ie 



(72) 



Thus the horizon itself possesses only axial symmetry. The deviation from spherical symmetry is small, though. For 
instance for the solution of Figs. 8 we find Le/Lp ~ 0.998. 



B. xh dependence of the solutions 

Let us now study the dependence of the static axially symmetric solutions on the isotropic black hole horizon radius 
xjj. To be specific, we consider black hole solutions with 7 = 1, n = 2 and fc = 1, the parameters also employed in the 
solution shown in Figs. 8. In Fig. 9 the energy density of the matter fields is shown for several black hole solutions, 
with values of the isotropic black hole horizon radius = 0.01, 0.1, 1 and 10, as well as for the globally regular 
solution. As noted above, for larger values of xn the global maximum of the energy density of the black hole solutions 
resides on the z-axis at the horizon, while a local maximum is located on the p-axis away from the horizon. With 
decreasing xn, the maximum on the p-axis away from the horizon increases and becomes the global maximum, while 
the maximum on the z-axis at the horizon diminishes. At the same time, a pronounced minimum develops on the 
p-axis at the horizon, aggravating the angle-dependence of the energy density at the horizon. 

With decreasing xh the energy density of the matter fields of the black hole solutions tends increasingly towards 
the energy density of the globally regular solution, which possesses a toroidal shape because of the strong global 
maximum on the p-axis |^,^^. However, the limit xh — > is not smooth, as already observed in section III for the 
static spherically symmetric solutions. The reason is, that the magnetic field of the black hole solutions is purely 
radial at the horizon, while the magnetic field of the globally regular solutions also has non- vanishing Bg at the origin. 
In the static axially symmetric solutions, both _B,. and Bg are angle-dependent. For the globally regular solutions the 
contributions from both Br and Bg precisely add to an angle- independent density at the origin. In contrast, the black 
hole solutions possess an angle-dependent density at the horizon. 

In Figs. 10-12 we show the functions of the black hole solutions with 7i = 2, fc = l, 7 = 1 and xh ~ 0.01, 0.1, 1 
and 10 for three angles. Figs. 10 show the metric functions. Figs. 11 the gauge field functions and Fig. 12 the dilaton 
function. For small values of xh, the angle-dependence of the metric functions is strongest around x — 1. With 
increasing xn, the angle-dependence decreases strongly, and the metric becomes increasingly spherical. At the same 
time, with increasing xh the matter fields become less important. This is seen for instance in Fig. 13, where we 
present relation (^), = po + 2TS' (with po = D/j), and observe that po ^ 2TS for small xn, whereas po ^ 2TS 
for large xh- Considering the gauge field functions, we recall, that for the globally regular solutions the functions H2 
and H4 have precisely k nodes, the function Hi has fc — 1 non-trivial nodes, and the function H3 has one non-trivial 
node ||]. For the black hole solutions we here observe accordingly, that Hi has no node and that H2 and H4 each 
have one node. However, we observe two non-trivial nodes for the function H3 of the black hole solutions of Fig. 11c, 
which may be due to the choice of gauge. We note, that the gauge field functions of Figs. 11 are obtained with gauge 
condition (^9|). With increasing xn the gauge field functions retain a considerable angle-dependence. They approach 
a limiting shape, shifting towards larger values of x. The dilaton function is slightly angle-dependent at the horizon, 
as seen in Fig. 12. Analogously to the metric functions, for small values of xn the angle-dependence of the dilaton 
function is strongest around x = 1. With increasing xn the angle-dependence decreases strongly, and at the same 
time the magnitude of the dilaton function diminishes strongly. 

In Fig. 14 we exhibit the Kretschmann scalar for these solutions. With decreasing xn the Kretschmann scalar of the 
black hole solutions tends to the Kretschmann scalar of the globally regular solution except close to the horizon, where 
it increases dramatically. For finite xn, however, the Kretschmann scalar remains finite at the horizon, indicating 
that the black hole solutions indeed possess a regular horizon, as required (see Appendix A.C). 



11 



In general the EYM solutions are very similar to the EYMD solutions. We therefore do not exhibit these here. For 
instance, the energy density of the matter fields and the metric functions / and m for the EYM black hole solutions 
with n = 2, k — 1 and xh = 0.02, 0.1, 0.5 and 1 are shown in Q]. 

Let us now inspect the metric functions at the horizon more closely. In Figs. 15a-c we exhibit the expansion 
coefficients f2{0), m2{9) and hiO) of the metric functions for the above set of EYMD solutions. With increasing xh 
the angle-dependence of the expansion coefficients first increases and then decreases again. The shape of the horizon 
changes with xh in a similar way. This is seen in Fig. 16, where the ratio of the circumference at the equator and 
the circumference at the poles, Le/Lp, is shown as a function of the mass for the EYMD (7 ~ 1) and EYM black 
hole solutions with n = 2 and fc = 1. The maximal deformation of the horizon occurs for an isotropic horizon radius 
of xu = 0.295 for 7=1 and x-a — 0.195 for 7 = 0. For 7 = 0, the maximal deformation of the horizon is greater 
than for 7=1. Fig. 17 shows the area of the horizon as a function of the mass for the corresponding solutions. For 
comparison, the area of a spherical horizon with circumference Lg is also shown, deviating only little from the area 
of the deformed horizon. 

Fig. 18 shows the inverse temperature as a function of the mass for the static axially symmetric EYMD (7 = 1) and 
EYM black hole solutions with n — 2 and k = 1. These curves are very similar to those of the corresponding static 
spherically symmetric black hole solutions ||2^. Table 1 presents the dimensionless mass, temperature, entropy and 
dilaton charge as well as the ratio L^/Lp of the above set of static axially symmetric black hole solutions of EYMD 
theory (7 = 1). 

C. Winding number dependence 

To illustrate the winding number dependence of the black hole solutions we show in Figs. 19 the solutions with 
n = 4, A:==l, 7=1 and xh = 0.01, 0.1, 1 and 10 for three angles. In Fig. 19a we see the energy density of the 
matter fields. We recall, that with increasing winding number n the maximum on the p-axis of the energy density 
of the globally regular solutions shifts outward and decreases in height As compared to the n = 2 black hole 
solutions, shown in Fig. 9, we here observe that with decreasing xh the globally regular solution is approached faster 
for the greater winding number n = 4, while with increasing xn the angle-dependence of the energy density of the 
black hole solutions remains stronger. Most strikingly, however, we observe, that the global maximum always resides 
on the /9-axis. The maximum on the z-axis remains a local one also for large xn. 

In Figs. 19b-d we show the metric function /, the gauge field function H2 and the dilaton function ip, respectively. 
For the globally regular solutions the angle-dependence of the metric and matter functions increases strongly with n 
and the location of the biggest angular splitting moves further outward. This dependence is refiected in the black hole 
solutions. In particular we observe, that with increasing n the angular dependence remains stronger for larger xh- 

Let us now turn to the shape of the horizon. As for n = 2, with increasing xu the ratio of the circumference at 
the equator and the circumference at the poles, L^/Lp, first decreases and then increases again. In Fig. 20 the ratio 
Le/Lp is shown as a function of the mass for the EYMD (7 — 1) and EYM black hole solutions with n = 4 and k — 1. 
The maximal deformation of the horizon occurs for an isotropic horizon radius of xn = 0.95 for 7 = 1 and xn = 0.65 
for 7 = 0. Again, for 7 = the maximal deformation of the horizon is greater than for 7=1, but for n = 4 it is 
smaller than for n = 2. Table 1 presents the ratio L^/ Lp as well as the dimensionless mass, temperature, entropy and 
dilaton charge for the above set of static axially symmetric black hole solutions of EYMD theory (7 = 1). 

D. Node number dependence 

To illustrate the node number dependence of the black hole solutions we show in Figs. 21 the solutions with n = 2, 
/c = 2, 7 = 1 and the same set of radii xjj as above. The energy density of the matter fields is shown in Fig. 21a. 
We recall that the maximum of the energy density of the globally regular solutions is located on the p-axis. With 
increasing node number k it shifts inward and increases strongly in height . For the black hole solutions we observe 
that with decreasing xn the globally regular solution is approached more slowly for greater node number k, while with 
increasing xn the angle-dependence of the energy density of the black hole solutions diminishes faster. In particular, 
for large xn the maximum on the z-axis becomes the global maximum, while the maximum on the p-axis disappears. 

In Figs. 21b-d we show the metric function /, the gauge field function H2 and the dilaton function if, respectively. 
We recall, that for the globally regular solutions the location of the biggest angular splitting moves inward with k 
while its size stays roughly constant for most functions M. For the black hole solutions we here observe, that as 
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compared to the k = 1 solutions, shown in Figs. 10-12, the angular dependence of the metric and dilaton functions 
diminishes faster with increasing xn. 

Table 1 again shows the dimensionless mass, temperature, entropy and dilaton charge as well as the ratio L^/Lp 
for the above set of static axially symmetric black hole solutions of EYMD theory (7 — 1). 



E. Limiting solutions 

For fixed n, 7 and xn and increasing k, the static axially symmetric black hole solutions form sequences, tending 
to limiting solutions. Whereas the solutions of the sequences are magnetically neutral, axially symmetric and non- 
abelian, the limiting solutions possess magnetic charge n and they are spherically symmetric and abelian. For finite 7 
the limiting solutions are EMD black hole solutions ||l^ , while for 7 = the limiting solutions are Reissner-Nordstr0m 
solutions jj,^ . The convergence of the global properties is seen for instance in Table 1 of . 

To illustrate the convergence of the metric and matter functions, we exhibit in Figs. 22a-c the EYMD solutions 
for n = 2, 7 = 1 and xh = 1 as an example. With increasing k, the metric functions converge rapidly to the metric 
functions of the limiting EMD solution. In Fig. 22a this is seen for the metric function /. With increasing k, the 
gauge field functions tend to their (vanishing) limiting functions in an exponentially increasing interval, but because 
of the boundary conditions the convergence is not uniform for H2 and H4. This is illustrated in Fig. 22b for the gauge 
field function H2. The dilaton function shown in Fig. 22c again converges rapidly and uniformly. 



V. CONCLUSIONS 

We have constructed numerically a new class of black hole solutions in E YM and EYMD theory jj] . These black 
hole solutions are asymptotically flat, static and possess a regular event horizon. However, they are not spherically 
symmetric but only axially symmetric with angle-dependent fields at the horizon. 

The event horizon of the static axially symmetric black hole solutions resides at a surface of constant isotropic 
radial coordinate, x = xn. However, the horizon is not spherical. Evaluating the circumference of the horizon along 
the equator, L^, and the circumference of the horizon along the poles, Lp, we observe that the ratio L^/ Lp is slightly 
smaller than one, i.e. the horizon is slightly elongated along the symmetry axis, the maximal elongation occurring for 
small values of xn. 

Like their globally regular counterparts, the static axially symmetric black hole solutions are characterized by two 
integers, the winding number n > I and the node number k of the purely magnetic gauge field. Whereas the energy 
density of the globally regular solutions has a torus- like shape, due to a strong peak on the p-axis away from the origin 
1^, the energy density of the black hole solutions has a more complicated shape, depending on the winding number 
n, the node number k and the horizon radius xjj. 

The static spherically symmetric EYM and EYMD black hole solutions are unstable |2^,|l^ , and there is all reason 
to believe, that the static axially symmetric black hole solutions are unstable, too. But we expect analogous black 
hole solutions in EYMH theory |^,^ and ES theory ||2^, corresponding to black holes inside axially symmetric 
multimonopolcs and multiskyrmions, respectively, and for n = 2 these static axially symmetric solutions should be 
stable |p|, ^8|j2S[| . In contrast, stable black hole solutions with higher magnetic charges (EYMH) or higher baryon 
numbers (ES) should not correspond to static axially symmetric solutions with n > 2. Instead these stable black hole 
solutions should exhibit only discrete crystal- like symmetries [|^,^,|9|. We expect analogous but unstable black hole 
solutions with crystal-like symmetries also in EYM and EYMD theory. 

We would like to thank the RRZN in Hannover for computing time. 
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VI. APPENDIX A 



A. Expansion at the horizon 



Here we present the expansion of the functions of the static axially symmetric black hole solutions at the horizon 
Xu in powers of S, 



6= 1 . 



(73) 



The expansion of the functions at the horizon can be obtained from the regularity conditions imposed on the Einstein 
equations and the matter field equations. 
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(74) 



The expansion coefficients /2, ^2, ^^ii, H20, H^q, H^q and ipo are functions of the variable 9. The expansion 
depends on the gauge condition imposed on the gauge field functions at the horizon. For the expansion above we 
employed the gauge condition (pj|), {deHi)\^^^^ — 0. 
The expansion of the functions at infinity is given in |9|. 



B. Surface gravity 

Here we show that the surface gravity Kgg of the static axially symmetric black hole solutions is indeed constant at 
the horizon. We employ the expansion of the metric functions of the static axially symmetric black hole solutions at 
the horizon xh in powers of S, eq. ([73|), 
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(75) 

Expanding the rd component of the Einstein equations at the horizon, we obtain a relation between the expansion 
coefficients f2{d) and m2{0), 

^ _ deni2 _ ^ 80/2 
m2 f2 

(76) 

Comparison with eqs. (p2h, (BSh then shows that the surface gravity is constant at the horizon. 



C. Kretschmann scalar 



Here we derive the expression for the Kretschmann scalar K, eq. ([G^), for the static axially symmetric black hole 
solutions and show that it is finite at the horizon. 

The nonvanishing components of the Ricmann tensor are 
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The Kretschmann scalar K = R^vXpR^^^'^ is then given by 
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In order to show that the Kretschmann scalar is finite at the horizon, we insert the expansion of the metric functions 
eq. ( [tbI ) into the Riemann tensor. In terms hke — --^) and {^j^ ~ the divergences cancel, because the 

lowest terms in the expansion are of the same order in r — rpj for all metric functions. The terms (2 + ^'-7^ — "f"^) and 
(2 + r^-f- — r^) are proportional to r — rn, and the term [r'^Li^ -|- 4r^ — r^^^ + 2) is finite at the horizon. The only 
term left, which could be divergent at the horizon, is (r^^^ — 2r^-^). However, this becomes _2»_(!IiM. _ 2^%^), 
which does not diverge because of relation ([76[). 
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VII. APPENDIX B 



Here we discuss our final choice of functions for the numerical integration. 

At the horizon the expansion coefficients f2{S), m2{0) and hiO) are of particular interest, because they enter into 
the expressions for the temperature, the area and the circumferences. In order to obtain these functions directly we 
have introduced the new functions f{x,9), fh{x,6) and l{x,6), 

f{x,e) =x^ f{x,e) , m{x,e) = x^ rfi{x,0) , l{x,0) = x^ l{x,e) , (77) 

where x = (1 — ^) is the compactified coordinate, eq. (|69|). The functions ./2(^^), ^2(0) and hi^) are then given by 

f2i0) = ^ /(O, e) , m2(9) = ^ m{0, 9) , UO) = \ [(0, 9) . (78) 

Xjj Xjj Xjj 

In the limit x ^ 00 the variable x approaches the value 1. Consequently, at infinity the boundary conditions for 
the new functions /, rh and I coincide with the boundary conditions of the functions /, m and Z, respectively. At the 
horizon the boundary conditions for the new functions can be obtained from the expansion of the metric functions. 
They are given by 

(/-^s/)L=o = 0' ('^ + 52m)U^o = 0, {Udd%^^^0. (79) 
To satisfy the regularity condition (^ exactly in the numerical calculations, we have introduced the function g{x, 9), 

fh(x,9) 

= (80) 
On the symmetry axis this function satisfies the boundary condition 

9\e=o = 1 (81) 

and at the horizon 

5s5U=o = . (82) 
The numerical calculations were mostly performed with the functions /(x, 9), g{x, 9) and l{x, 9). 
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EYMD(7 = 1) 




k=l, n=2 


k=l, n=4 


k=2, n=2 




M 


0.01 


0.9766 


1.6227 


1.2762 


0.1 


1.1183 


1.7700 


1.4081 


1. 


2.5649 


3.2193 


2.5649 


10 


20.079 


20.244 


20.079 




T 


0.01 


0.583 


0.989 


0.230 


0.1 


0.727 


0.076 


0.039 


1. 


0.016 


0.013 


0.016 


10. 


0.002 


0.002 


0.002 




S 


0.01 


0.014 


0.012 


0.033 


0.1 


1.167 


1.116 


2.243 


1. 


62.47 


71.25 


62.47 


10. 


5045. 


5081. 


5045. 




D 


0.01 


0.9609 


1.6063 


1.2619 


0.1 


0.9487 


1.6011 


1.2319 


1. 


0.6022 


1.3188 


0.6023 


10. 


0.0836 


0.2682 


0.0836 




Le. / 


0.01 


0.9999 


1.0000 


0.9993 


0.1 


0.9974 


0.9997 


0.9910 


1. 


0.9980 


0.9968 


0.9980 


10. 


1.0000 


0.9998 


1.0000 



Table 1 

The dimensionless mass /x, the temperature T, the entropy S, the dilaton charge D and the ratio L^jLp of the EYMD 
black hole solutions with node number k= \ and winding number n = 2, A; = 1 and n = 4, as well as A; = 2 and n = 2 
are shown for the values of the horizon radius xh = 0.01, 0.1, 1 and 10. 
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n=1, y=1, x^=0.02 
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Fig. 1 

The coordinate transformation between the isotropic coordinate x and the Schwarzschild-like coordinate x is shown 

for the static spherically symmetric solutions in = 1) of EYMD theory with 7=1 and = 1 — 3. Also shown is the 
coordinate transformation for the limiting EMD solution. 
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Fig. 2a 

The metric function / is shown for the static sphericaUy symmetric solutions (n = 1) of EYMD theory with 7=1 
and A; = 1 — 3. Also shown is the metric function of the limiting EMD solution. 
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Fig. 2b 

Same as Fig. 2a for the metric function m. 
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Fig. 3a 

Same as Fig. 2a for the gauge field function w, for /c = 1 — 4. 
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Fig. 3b 

Same as Fig. 2a for the dilaton function ip. 
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Fig. 4 

The coordinate transformation between the isotropic coordinate x and the Schwarzschild-like coordinate x is shown 
for the static spherically symmetric solutions (n = 1) of EYM theory with fc = 1 — 3. Also shown is the coordinate 
transformation for the limiting RN solution. 
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Fig. 5a 

The metric function / is shown for the static spherically symmetric solutions (n = 1) of EYM theory with fc = 1 — 2. 
Also shown is the metric function of the limiting RN solution. 
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Fig. 5b 

Same as Fig. 5a for the metric function m for /c = 1 — 3 . 
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Fig. 6a 

The energy density e = —Tq is shown for the static sphericaUy symmetric black hole solutions (n = 1) of EYMD 
theory with 7 = 1 and k = 1 and the horizon radii xu = 1, 0.1, 0.01 and 0.001, as well as for the corresponding 
globally regular solution. 




Fig. 6b 

Same as Fig. 6a for k = 2. 
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Fig. 7 

The Kretschmann scalar is shown for the static spherically symmetric black hole solutions (n = 1) of EYMD theory 
with 7=1 and k = 1 and the horizon radii Xu = 1, 0.1, 0.01 and 0.001, as well as for the corresponding globally 
regular solution. 
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Fig. 8a 

The energy density e = — Tq is shown for the static axially symmetric black hole solution of EYMD theory with 7 = 1, 
n = 2 and k = 1 and the horizon radius a;H = 1 in a 3-dimensional plot and a contour plot with axes p and z. 
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Figs. 8b-e 

Surfaces of constant energy density e = — Tq are shown for the solution of Fig. 8a. 
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Fig. 9 

The energy density e = — Tq is shown as a function of the dimensionless coordinate x for the angles ^ = 0, ^ = 7r/4 
and = 7r/2 for the EYMD black hole solutions with 7=1, winding number n = 2, node number fc = 1 and horizon 
radii = 0.01, xh = 0.1, xh = 1, and a;H = 10, as well as for the corresponding globally regular solution. 
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Fig. 10a Same as Fig. 9 for the metric function /. 
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Fig. 10b Same as Fig. 9 for the metric function m. 
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Fig. 10c Same as Fig. 9 for the metric function I. 
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Fig. lid 



Fig. 11a Same as Fig. 9 for the gauge field function Hi 

Fig. lib Same as Fig. 9 for the gauge field function H2 for horizon radii xu = 0.1, xu = 1, and xh = 10. 
Fig. 11c Same as Fig. 9 for the gauge field function H3. 
Fig. lid Same as Fig. lib for the gauge field function H4. 
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Fig. 12 

Same as Fig. 9 for the dilaton function y>. 
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Fig. 13 

The dimensionless mass /x is shown as a function of the horizon radius ojh for the EYMD black hole solutions with 
7 = 1, winding number n = 2, node number fc = 1. Also shown are 2TS and -D/7. 
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Fig. 14 

The Kretschmann scalar is shown as a function of the dimensionless coordinate x for the angles 6 = 0, ^ = 7r/4 and 
= 7r/2 for the EYMD black hole solutions with 7 = 1, winding number n = 2, node number k = 1 and the horizon 
radii = 1, 0.1, and 0.01, as well as for the corresponding globally regular solution. 
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Fig. 15a 

The normalized expansion coefficient /2 of the metric function / is shown as a function of the angle 9 for the EYMD 
black hole solutions with 7=1, winding number n = 2, node number k = 1 and the horizon radii xa = 10, 1, 0.1, 
and 0.01. 
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Fig. 15b Same as Fig. 15a for the expansion coefficient 1112 of the metric function m. 
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Fig. 15c Same as Fig. 15a for the expansion coefficient li of the metric function I. 
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Fig. 16 

The ratio Lg/Lp of the circumference of the horizon along the equator Lg to the circumference of the horizon along 

the poles Lp is shown as a function of the dimcnsionlcss mass /i for the black hole solutions with winding number 
n = 2 and node number A; = 1 of EYMD theory with 7 = 1 and EYM theory. 
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Fig. 18 

Same as Fig. 16 for the inverse temperature (3 = T' 



37 



Fig. 19a 





log,c(x) 



0.0 



2 
Log,„(x) 



Fig. 19b 



ic=1 , n=4, y = 1 
x„=0.1 ^ x„=1 

... \ x„=10 



Fig. 19c 



2 
log,„( x) 



e=7T./4.. - 
e=Tj/z_ 





-. 1 

-.2 



k=1 , n=4, y = 1 



• x„=10 — 

"7 


e= 

e=TT/4 " 


x„= 1 y/'l 

1 


e=Tj/z_ _ 


x„=0.1 ,,.^>flil 




- reg ^_^y'~ \ 





2 
log,„( x) 



Fig. 19d 



Fig. 19a 

The energy density e = — Tq is shown as a function of the dimensionless coordinate x for the angles 6 = 0, 6 = tt/A 
and 6 = w/2 for the EYMD black hole solutions with 7 = 1, winding number n = 4, node number k = 1 and horizon 
radii xn = 0.1, xn = 1, and xn = 10, as well as for the corresponding globally regular solution. 



Fig. 19b 

Same as Fig. 19a for the metric function /. 



Fig. 19c 

Same as Fig. 19a for the gauge field function H2. 



Fig. 19d 

Same as Fig. 19a for the dilaton function (fi. 
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Fig. 20 

The ratio Lf./Lp of the circumference of the horizon along the equator Lg to the circumference of the horizon along 

the poles Lp is shown as a function of the dimcnsionlcss mass fi for the black hole solutions with winding number 
n = 4 and node number A; = 1 of EYMD theory with 7 = 1 and EYM theory. 
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Fig. 21a 

The energy density e = — Tq is shown as a function of the dimensionless coordinate x for the angles 6 ^ 0, 9 = tt/A 
and 6 = w/2 for the EYMD black hole solutions with 7 = 1, winding number n = 2, node number k = 2 and horizon 
radii xu = 0.01, xn = 0.1, xn = 1, and xa = 10, as well as for the corresponding globally regular solution. 



Fig. 21b 

Same as Fig. 21a for the metric function /. 



Fig. 21c 

Same as Fig. 21a for the gauge field function H2. 



Fig. 21d 

Same as Fig. 21a for the dilaton function (fi. 
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Fig. 22a The difference of the metric function and the metric function /oo of the hmiting sohition is shown as a 
function of the dimensionless coordinate x for the angles — 9 = -k /A and 9 ^ 11/2 for the EYMD black hole 
solutions with 7 = 1, horizon radius xyl — 1, winding number n = 2, and node numbers & = 1 — 3. 
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Fig. 22b The gauge field function H2 is shown as a function of the dimensionless coordinate x for the angles 9 = Q, 
9 = 7r/4 and = 7r/2 for the EYMD black hole solutions with 7=1, horizon radius xh = 1, winding number n = 2, 
and node numbers fc = 1 — 3. Also shown is the gauge field function of the limiting EMD solution. 




Fig. 22c Same as Fig. 22b for the dilaton function {p. 
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